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Abstract 

(N 

^ I Regge theory provides a very simple and economical description of data for (i) the proton 

'nI" . structure function with x < 0.07 and all available values, (ii) the charm structure 

function, and (iii) jp ^ J/ipp. The data are all in agreement with the assumption that 
\^ ' there is a second pomeron, with intercept about 1.4 . They suggest also that the contribution 

CD ■ from the soft pomeron is higher twist. 

This means that there is an urgent need to make perturbative evolution compatible with 
Regge theory at small and not-so-small x, and to reassess the magnitude of higher- twist 
contributions at quite small x. 
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^ . Regge theory is one of the great truths of particle physics^^]. It is supposed to be applicable if W"^ 
d \ is much greater than all the other variables. Thus, we expect it to be valid when x is small enough, 
even for large values of . Fits to large-Q^ data have traditionally concentrated on perturbative 
evolution; in this paper we report fits that rather emphasise the Regge behaviour. Regge behaviour 
is not a substitute for perturbative evolution, but a constraint on it, and an important task will be 
to understand better how the two coexist. Somewhat surprisingly, we find that Regge behaviour is 
compatible with the data even for values of x up to 0.07 or higher. 

Regge theory involves pomerons. Until HERA measurements of the proton structure functiont^l at 
very small-x, and of J/V' photoproductionl^l , it seemed consistent to suppose^ that a single nonper- 
turbative pomeron with intercept close to 1.08 describes the whole of diffractive physics, including!"'] 
the NMC and E665 structure function data at fairly small x. The more violent behaviour observed 
at HERA calls rather for an intercept closer to 1.5 than 1.0, and there has been much discussion of 
what is responsible for this. 

One viewl^] is that the pomeron intercept varies with , but we find this difficult to believe. Regge 
theory relates'^] the Vl^-dependence of a process to the positions of singularities in the complex angular 
momentum plane, and it does not naturally accommodate the notion that these positions vary with 
. The relative weights of the contributions from these singularities can, however, vary with Q^, 
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and so the effective power resulting from combining them can be Q^-dependcnt. We find less than 
persuasive the hypothesis'^] that in soft processes there is a significant amount of shadowing, so that 
the observed effective intercept of 1.08 actually is the result of the pomeron having a significantly 
larger intercept, and that the rapid rises observed at HERA are merely the result of the disappearance 
of the shadowing. If there is a significant amount of shadowing, the apparent intercept should be 
process-dependent, but all soft processes are foundlS], to a remarkable accuracy, to have the same 
value. 

In this paper, therefore, we explore the notion that there are two pomerons. It is a good working 
principle always first to try the simplest assumptions, so we assume that each is a simple pole in 
the complex angular momentum plane, making the contribution from each a simple power of W, 
multiplied by an unknown function of Q^. Thus 

i 

Bearing in mind that the cross section for the absorption of real photons is 



iir aEM rn 



(2) 

Q2=0 



we require each F2{x, Q^) to vanish linearly with at fixed W, so that presumably fi{Q^) vanishes 
like (Q'^y~^'^\ Beyond the fact that, presumably, they are smooth functions of Q^, this is the only 
information that we have about the form of the fi{Q^)- In the past 1^1 , we used the simplest possible 
forms: 

fiiQ') = Ai(^^y^^^ (3) 



and just two terms, with powers 



ei = 0.0808 (soft pomeron exchange) 

€2 = -0.4525 (/, a exchange) (4) 

Now, we retain the same two powers ei and 62, and we add a third term (i = 0), which we call 
hard-pomeron exchange. 

As we have said, Regge theory is supposed to be applicable if W'^ is much greater than all the other 
variables. Thus, we expect it to be valid when x is small enough, whatever the value of Q^. Just what 
is meant by small enough is open to debate: we shall assume, as we have beforel^l , that it means 



X < 0.07 



and, in the case of the real-photon data, ^/s > 6 GeV, though it is worth remembering that these 
choices are somewhat arbitrary. 

Because Regge theory gives us no information about the dependence of the various contributions, 
we should not expect our simple guess (3) for the coefficient functions /i(Q^) to work for arbitrarily 
large Q^. So initially, again as beforel^l, we used it to extract a preliminary value for the power €0 by 
applying it to all data for which 

< 10 GeV^ 

including = 0. Wc then abandoned the simple form (3), and used this preliminary value for eo to 
fit at each value of Q^, and so found how the data require the coefficient functions fi{Q'^) to behave at 
all available values of Q^. This suggested how we should parametrise the functions fi{Q^), and with 
these parametrisations we then began the fitting exercise afresh. 
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Figure 1: Best fits to the coefHcient functions /o(Q^) and /i(Q^), at each for which there are data 
with X < 0.07, together with curves corresponding to the fit (6) and (7). 

The simple form (3) gave = 1.075 per data point*, for the 280 data points with < 10, with 
parameter values 



Ao = 0.0358 
Ai = 0.274 
A2 = 0.149 



eo = 0.435 



ao = 5.41 
ai = 0.504 
02 = 0.0196 



(5) 



We have given each parameter to high accuracy because a small change in any one parameter has 
a significant effect on the of the fit. This does not mean that the parameters are determined to 
that accuracy, because the minimum is not very sharp and a change in any one parameter can be 
compensated by changes in all the others in such a way as not to have much effect on the x^- 

For larger Q"^, we find that the simple dependences make the coefficient function of the x~'^° term 
too small, and that of x~^^ too large. For each value of for which we have data with x < 0.07, we 
made a best fit using the form (1) keeping the same simple shape (3) for the /, a contribution and with 
the values of all the quantities except /o(Q^) and fi{Q^) taken over from the fit to the < 10 data. 
The resulting values for fo{Q^) and fi{Q^) are shown in figure 1, where the error bars are produced 
from MINUIT. Obviously one can fit these with a variety of functions. We have tried several. Our 
lowest overall corresponds to the choices 



+ ao 



/i(Q') = ^i 



l+eo 



l + Xlog(l + ^^ 



h{Q^) = A 



l + ei 



1 



1 + 

^2 \ 1+^2 



+ 02 



(6) 



We vary all the parameters except ei and €2. This 10-parameter fit gave a minimum %^ = 1.016 per 
data point for 539 data points with the parameter values 



* In our calculations of we have folded the systematic and statistical errors in quadrature. 
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eo = 0.418 

Ao = 0.0410 ao = 7.13 

X = 0.485 Ql = 10.6 

Ai = 0.387 ai = 0.684 Qj = 48.0 

A2 = 0.0504 aa = 0.00291 (7) 

The same comments as before must be made about the accuracies of these parameters. Figures 2, 3 
and 4 show the fits they give to the data. 



We have a number of comments to make: 

1. Our fits show that, at small enough x, the data may be fitted by a combination of three fixed 
powers of x. This is in accord with the Regge-theory expectation if, at least to a good approximation, 
the singularities in the complex angular momentum plane are simple poles. We have used as input the 
same powers, ei = 0.0808 and €2 = —0.4525, for soft pomeron and /, a exchange as in our previous 
fitsl^l, though these values have recently been questioned by Cudell et al^^^l. Instead of (4), they 
advocate ei = 0.094 and €2 = —0.33. With these values, and the form (6) for the coefficient functions 
fi{Q^), our per data point would reduce to 0.98, with eo = 0.410, though whether this reduction 
is significant is doubtful, since our is already so close to 1. 

2. There may be a significant error on our output hard-pomeron power eo =0.418, perhaps ±0.05. For 
example, if one were to decide that the value 02 = (53.9 MeV)^ in (7) is too small to be physically 
reasonable, then increasing it to the square of the pion mass increases the per data point to 1.032, 
and so still gives a good fit, but it reduces eo to 0.388. Alternatively, if we prefer to have the coefficient 
fu nction f i{Q^) for the soft pomeron go to zero at large Q"^ like l/Q^ instead oil/Q, so that we replace 
^JQ^IQl in (6) with Q^/Qf , the per data point becomes 1.018 with eo = 0.462. 

3. This makes it clear that the exact forms we have used for the dependences of the three 
contributions should not be taken too seriously: Regge theory has nothing to say about them and we 
have chosen the simplest forms that work. Nevertheless, we can be reasonably confident that, as is 
seen in figure 1, the soft-pomeron term has a coefficient that rises to a maximum somewhere between 

= 5 and 10, and it is very likely that it goes to zero at high Q^. The hard and soft pomeron 
contributions to arc shown in figure 5 at two values of they move in opposite directions as 
increases, so that while the soft pomeron is important at = 5 it becomes less so at higher . The 
conclusion that the soft pomeron is higher twist is one that we did not expect and implies that all fits 
based on conventional evolution that have been made so far need to be reassessed. It is reminiscent of 
a suggestion of McDermott and Forshaw^^^] , who indeed have a ncxt-to-leading contribution that goes 
to zero as l/Q at high . Their leading term rises like Q at high the data are not inconsistent 
with this: we find it gives a per data point of 1.072, with eo = 0.385. But our logarithmic form, 
which presumably corresponds more to DGLAP evolution, gives a better fit though we emphasise 
that Regge theory requires that the perturbative evolution should be applied in such a way that, at 
small it does not cause the power eo to change with Q'^ . 

4. We have said that our choice x < 0.07 for the range of x in which we fit the data is arbitrary. Indeed, 
one might be surprised that the Regge region extends to such a large value of x. If we use instead 
only data for which x < 0.01, with the functional forms (6) for the coefficient functions fi{Q'^) and 
the parameter values (7), the x^ improves from 1.016 per data point to 0.95 (with 357 data points). 

5. There will be some who will object to fixed powers of x because of worries about unitarity. They 
will be wrong to do so, for several reasons. Firstly, we have shown that the fixed powers fit the data, 
and the data certainly respect whatever constraints unitarity may impose. Secondly, unitarity has 
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Figure 2: Fits to real-photon and small-Q^ data. In the small-Q^ plots the vertical axes are i^'Ji and iFj, so that the larger x-values 
are at the top. 
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Figure 3: Fits to F2 data for between 0.11 and 10. Tiie vertical scales are F^IQ'^, so the smaller-Q^ data are at the top 
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Figure 4: fits to large-Q^ data. The vertical scales are F^jQ'^, so the smaller values are at the top 




Figure 5: Relative magnitudes of hard and soft pomeron contributions at = 5 and 1000 



nothing useful to say about structure functions anyway: the derivation of the unitarity constraints on a 
purely hadronic process such as pp scattering relies on the fact that a pp intermediate state appears in 
the unitarity equation; but 7*p scattering is very different, because there is no 7*p intermediate state. 
In any case, even in pp scattering, data at present energies are not seriously affected by unitarity 
constraints. The fact that a constant power fits the data well^ indicates that, so far, shadowing 
corrections to the t = amplitude are small, though they will become more important at higher 
energies. Similarly, at extremely small x the fixed-power behaviour will be moderated by shadowing 
suppression, but there is no reason to believe that such small values have yet been achieved. 

6. The accurate data for the low-energy real-photon cross-sections are an important constraint. If we 
make a best fit to the virtual-photon data alone, its extrapolation comes nowhere near the low-energy 
real-photon data. As can be seen in figure 2, our new fit extrapolates the small-Q^ data to give a 
larger real-photon cross-section at HERA energies than in our old fitsW. The ZEUS collaboration 
has reached a similar conclusiont^^l , but it must be emphasised that it depends crucially on assuming 
a simple form for the dependence. There is no reason to rely on this, and by adopting more 
complicated extrapolations one can arrive at lower valuesl^"^]. We could even decouple the hard 
pomeron completely at = 0. Nevertheless, it is interesting to note that preliminary results from 
a new ZEUS measurement of the real-photon total cross-section agree rather well with our curve^^^l , 
which is 

i 

= 0.283 (2z/)°-4is + 65.4 (2z/)°-°8°8 + 138 (2i/)-°-4^25 



with 2i/ in GeV^ and the coefficients in /xb, where 2u = — m^. 

7. We have made the assumption that the Q"^ dependences of the contributions to F2 from different 
flavours are the same, or at least that when the contributions are added together their sum has simple 
behaviour. This is not obviously true, and it is interesting to subject the recent accurate preliminary 
ZEUS datal^^] for the charm content of F2 to our parametrisation. We find that they indicate that 
both the soft-pomeron and the /, a exchanges do not couple. The latter is readily understandable 
because of Zweig's rule, but the decoupling of the soft pomeron comes as something of a surprise. We 
impose 

^2' = ( j (1 + log(QVQc)) (10) 

with eo = 0.418 as before. The data are not sufficient to determine the remaining four parameters well, 
so we choose to impose a constraint that the coupling of the hard pomeron is flavour-blind at large 



8 



1e+09 
1e+08 
1e+07 
1e+06 
100000 
10000 
1000 
100 
10 
1 

0.1 



ZEUS CHARM (PRELIMINARY): Q'^2=2.5, 7, 11, 18 

30,60,130 



0.0001 



i 



0.001 



0.01 



1e+07 



1e+08 



100000 



o 



10000 



1000 



HI CHARM: Q"2=7, 12, 25, 45 



100 
0.0001 



0.001 

X 



0.01 



10 



0.1 



0.01 



0.001 



EMC CHARM: 0^^2=1 .33, 2.4, 4.2, 7.4, 13.3, 23.5 




' :::f r 







0.01 



0.1 



Figure 6: Fits to data for F^'^. The vertical scales are Q^Jb for the ZEUS and HI data, and Q'''F2 for EMC, so the smaller values 
are at the bottom for each plot. For the comparison with EMC data, we have scaled the fit by a factor 0.5. 
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Figure 7: Conjectured soft and hard pomeron trajectories, with 2++ gluebah candidatesl^^^ [^^1 

enough Q^. Assuming that so far is not large enough for there to be a significant contribution from 
b quarks, this implies that the coefficient of the logarithm in the hard-pomeron coefficient function for 
-^2 is |/(| + ^ + | + |) times that for the complete F2 : 

AcXc = OAAoX (11a) 

Then, the 3-parameter fit to the ZEUS data yields 

Ac = 0.184 ac = 8.45 Xq = 0.431 = 212 (116) 

and is shown in figure 6. This figure also includes comparison of the fit to the earlier HI and EMC 
datal^^l. In the case of the EMC data, we have had to renormalise our fit with a constant factor 0.5 to 
get good agreement. We cannot understand why this should be, but the latest MRS 

fittie] similarly is 

twice as large as the low-x EMC charm data. Ideally, one would first make a fit to the charm-quark 

data, then subtract it from the complete F2 and make another fit to just the light-quark contribution, 
but this must wait until the data improve. Eventually, the charm data will be important for reducing 
the error on the value of the hard-pomeron power eo: taken literally, the present charm data suggest 
that it should be towards the upper end of the range 0.42 it 0.05 that we have quoted above. 

8. We have found that we need a "hard-pomeron" term, with Regge intercept about 1.4. It is perhaps 
ironic that we reach this conclusion just as the hard BFKL pomeron seems to be in retreat [^^1. This 
does not rule out the possibility that our hard pomeron is perturbative. After all — as we have 
consistently maintainedl-^^l'I-^-^] — in the BFKL context asymptopia is very far away from present 
physics. Whatever the explanation of the hard pomeron, it is interesting that the mass scale Qq 
that determines how rapidly its contribution to F2 rises with is considerably larger than the 
corresponding soft-pomeron scale oi. Nevertheless, it may well be that even the hard pomeron is 
nonperturbative and that it is just a glueball trajectory. There is already some experimental support 
for the notion that this is true for the soft pomeron: the 2"'""'" glueball candidatet^^l at 1926 MeV is 
at just the right mass to lie on a straight trajectory of intercept close to 1.08 and slope 0.25 GeV~^. 
Another 2++ glueball candidate has just been announced , with a mass 2350 MeV. Let us suppose 
that this lies on the hard pomeron trajectory, and that this is also straight. The hard and soft 
pomeron trajectories will apparently then intersect in the complex angular momentum plane; however 
it is likely* that mixing will rather cause them to avoid each other in the way shown in figure 7. This 
picture gives the hard-pomeron trajectory a slope close to 0.1 — though we caution that figure 7 
should not be taken too seriously. 

9. Although data for quasi-elastic p photoproduction are well described by soft pomeron exchange, 
this is not the case for J/ip photoproduction. Not only is there a much more rapid rise with increasing 

* We are grateful to Peter Collins for a correspondence about this. 
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Figure 8: cr(7p — > J/ipp) in microbarns, together with simple fit 

energyl^^l , but also there is much less shrinkage of the forward peakl22]. This second feature suggests 
that the slope of the hard pomeron's trajectory may be small, at least near the forward direction, as 
we have indeed shown in figure 7. Assuming this, the contribution to the amplitude from the exchange 
of the two pomerons would be of the form 

1 1 
T{s,t) =i^Bi{t)s'^^''> e-2^'^«'(*) 

i=0 

eo(t) = eo + 0.1t ei(t) = ei + 0.25i (12a) 

The functions Bi{t) include products of the form factors that couple the pomerons to the proton and 
the jJ/ip vertex. We model them with the functions 

Bi{t) = B0)e^* (126) 

though their exact shape is not critical for our analysis - it has just a small effect on the energy 
dependence. The best fit to the existing data has 

5o(0)/5i(0) = 0.05 (12c) 

and is shown in figure 8. It is interesting that HI has new preliminary data points at higher energiesl^^l 
which agree well with our fit and confirm that the curve really should be concave upwards. 

10. An obvious question is whether the hard pomeron contributes also to purely soft processes. There 
is perhaps a hint of a small contribution to exclusive p photoproduction, in that the t-slope seems 
to be slightly smaller than what purely soft-pomeron exchange would give[23]. The contribution of 
the soft pomeron to the pp and pp total cross-sections was already fixed rather accurately I^^l by the 
data up to ISR energies, but there is room for some extra contribution in the higher-energy data. 
With some small downwards adjustment of the soft-pomeron term, one can tolerate a hard-pomeron 
contribution of perhaps 

0.01 (13) 

This would give slightly better agreement with the UA4 measurement than we had before^^^] but 
would still favour the E710 measurement over that of CDF[261. 
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11. The coefficient of the hard-pomeron term (13) in the pp total cross section is more than 2000 
times smaller than for the soft-pomeron term. However, if it turns out that (9) is the correct ^ 
extrapolation, then in the 7p total cross section the ratio is an order of magnitude larger. This would 
support the viewpoint of those^^^l who have maintained that the jp total cross section contains a 
significant piece that is different in character from purely hadronic cross sections. This additional 
piece cannot factorise in the Regge sense, because it would give a 77 total cross section much larger 
than has been measured l^^l . However, the preliminary OPAL and L3 data do hold the prospect that 
the hard pomeron may be visible in 77 collisions. 

12. The main message of this paper has been that Regge theory is applicable even at very large Q^, if 
only X is small enough. In this, our work does not relate immediately to any of the many other fits in 
the literature: we insist that the hard pomeron is present already at small and that evolution 
keeps both it and the other Regge powers intact. Kerley and Shawl^^l also have two pomerons, though 
the details of their work are very different. Our approach is in the same spirit as that of Haidtl^^l, 
though the details are again very different. The work of Nikolaev, Zakharov and Zollcrl"^*^' is closest to 
ours, though they do not extend their fit below = 1.5 and therefore have rather different powers of x 
from ours. Like us, they find that the coefficient functions of the nonleading x powers peak somewhere 
around = 10, though their behaviour at large is very different from ours. Indeed, we suggest 
that the soft pomeron contribution, although large for as large as 10 or more, is probably higher 
twist. 
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